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We reconsider the determinations of the strange quark mass m s from e + e _ into hadrons data using a new 
combination of FESR and revisiting the existing r-like sum rules by including non-resonant contributions to the 
spectral functions. To order and including the tachyonic gluon mass A 2 contribution, which phenomenologically 
parametrizes the UV renormalon effect into the PT series, we obtain the invariant mass rh s = (119 ± 17) MeV 
leading to: m s (2 GeV) = (104±15) MeV. Combining this value with the recent and independent phenomcnological 
determinations from some other channels, to order aj? and including A 2 , we deduce the weighted average: m s (2 
GeV)=(96.1 ± 4.8) MeV . The positivity of the spectral functions in the (pseudo)scalar [resp. vector] channels 
leads to the lower [resp. upper] bounds of m s (2 GeV): (71 ± 4) MeV < m s (2 GeV) < (151 ± 14) MeV, to 
order af. Using the ChPT mass ratio = 2m s /(m u + rrid) = 24.2 ± 1.5, and the average value of m s , we 
deduce: (m„ + m d )(2 GeV)=(7.9 ± 0.6) MeV, consistent with the pion sum rule result, which, combined with the 
ChPT value for m u /m d , gives: m d {2 GeV)=(5.1 ± 0.4) MeV and m u {2 GeV)=(2.8 ± 0.2) MeV. Finally, using 
(rn u +rri d ) from the pion sum rule and the average value of m s (without the pion sum rule), the method gives: 
r,3 = 23.5 ± 5.8 in perfect agreement with the ChPT ratio, indicating the self-consistency of the sum rule results. 
Using the value: mi,(mf,) = (4.23 ± 0.06) GeV, we also obtain the model-building useful scale-independent mass 
jj] ■ ratio: mf,/m s = 50 ± 3. 
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1 Introduction 



The determination of the strange quark mass is of prime importance for low-energy phenomenology, for CP- 
violation and for beyond standard model-buildings. Since the advent of QCD, where a precise meaning for the 
definition of the running quark masses within the M^-scheme £Q has been provided, a large number of efforts 
have been devoted to the determinations of the strange quark mass 1 using QCD spectral sum rules (QSSR) 2 a 
la SVZ 0, in the pseudoscalar 013110 EI GUI, the scalar El 112 , the e+e" [13 d U3 [H H3 H channels, 
tau-decay data 1191 f50] and lattice simulations [21U22II23II24| . while some bounds have been also derived from 
the positivity of the spectral functions [2H1 El ID IE! and from the extraction of the quark condensate [23 13 HI ■ 
In the following, we reconsider the determinations of m s from e + e~ into hadrons data by using new combinations 
of FESR ZU] and by revisiting the analysis done in El- I n so doing, we take into account more carefully the 
small non-resonant contributions into the spectral functions, though negligible in the individual sum rules become 
important in the combinations sensitive to leading order to m s . We also present a new combination of sum rule 
used in |2U| that we confront with previous sum rules presented in |15j . We conclude the paper by a comparison 
of recent different determinations of m s from QCD spectral sum rules from which we extract the average. This 
average being confronted to lattice calculations. 



2 Normalizations and notations 

We shall be concerned with the transverse two-point correlator: 

rC(<z) = * I d A x e^(0\TJ^(x) (J^O)) 1 10) = -(g^q 2 - q^)Tl ab (q 2 



J» EM = Vg{x) + ±=Vg{x) , (2) 



built from the SU(3) component of the local electromagnetic current: 

1 

7! 



where: 



V^{x)^\jU{x)\ a r^(x) ; (3) 



X a are the diagonal flavour SU(3) matrices: 



i / 1 \ . fT ' 1 



^ = \l2\ 1 ' As = V6 1 ' 1 (4) 



/ ' " \ -2 

acting on the basis defined by the up, down and strange quarks: 

/ u(x) \ 

i>(x) = d{x) . (5) 
V *(*) I 

In terms of the diagonal quark correlator: 

11% (q) i J d*x<*>*{0\TJ!f(x) (J;(0)) f |0> - -(g^q 2 - q^U^q 2 ) j = u,d,s, (6) 
where Jj(x) = tpjj^ipj, the previous SU(3) flavour components of the electromagnetic correlator read 3 : 
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(7) 



1 For reviews, see e.g. pi l3l 111 |5l . 

2 For a review see e.g. 151 . 

3 We shall follow the normalization used in SN 1151 . 



Therefore, the e + e — » hadrons total cross-section reads: 

4n 2 a o 1 



a(e + e — * hadrons) M 



d,s 



1 2 

<j Imn 33 (s) + -Imn 88 (s) + - 7 =Imn 38 (s) 



(8) 



In a narrow- width approximation (NWA), the resonance H contributions to the spectral functions can be intro- 
duced through: 



(oiki#; 



where the coupling ^Ha is related to the meson leptonic width as: 
which is itself related to the total cross-section: 



2 2 M H 
3" 



<j(e + e- -» H) = 12 ^ H ^ e+e ~ 5(s - M 2 ) 



(9) 



(10) 



(11) 



3 QCD corrections and RGI parameters 



In order to account for the radiative corrections, one introduces the expressions of the running coupling and 
masses. 

• To three- loop accuracy, the running coupling can be parametrized as |27ll5]: 



o.M = ai° ) <jl-a(°)|loglog A . J 



v 



(„,,) 3 [ | log Mo4-|i„ s i„4-| + |] + «}, 



with: 



,(o) 



-A log (v/A) 

and Pi are the 0(a\) coefficients of the /3-function in the A/S-scheme, which read for three flavours [Jp: 

Pi = -9/2, = -8, 3 = -20.1198. 
• The expression of the running quark mass in terms of the invariant mass is ^ [S]: 



(12) 
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a 2 s (v) + 1.95168af(i/) >, 



(15) 



where 7, are the 0(a*) coefficients of the quark-mass anomalous dimension, which read for three flavours |S]: 

7l =2, 72 = 91/12, 73 = 24.8404. (16) 
• The perturbative expression of the correlator reads, in terms of the running coupling evaluated at Q 2 = v 2 5 : 



l + U s = 



_ a s (Q 2 ) 



1.6398a: - 



10.2839- 



0l\ (* 



(17) 



where the last extra term in the c? s coefficient compared with the expression of the spectral function Im II comes 
from the analytic continuation. 

• The D = 2 contribution reads to order af, in terms of the running mass and by including the tachyonic gluon 
mass A 2 term [23 EH3 U EH1 : 



Q 2 Tl£= 2 \Q 2 ) 



1 

'4^2 



1.05a s A 2 + 6to 2 (1 + 2.6667a s + 24.1415a 2 + 250.4705a^ + C(a*)) 



(18) 



The coefficient of the a 4 term like the ones of all unknown higher order terms will be mimiced by the A 2 -term 
[HI 113 EH El EH] present in the D = 2 and D = 4 contributions. The presence of A 2 in the OPE helps in 
resolving the old puzzle of hierarchy scale [31] encountered in the sum rules analysis of the pion [SJ 0] E] and 
gluonia |35| channels. A 2 also improves the determination of a s and m s from r decay data |33M20| . In the present 
paper, the series converge slowly in the region where the analysis is performed. However, we expect that this 
slow convergence will not ruin the result, though introducing a large error, as each corrections are individually 
smaller than the lowest order term, while the size of the A 2 contribution (see Eq. (|19|l below) introduced to 
mimic the resummation of the unknown higher order terms remains a correction of the lowest order one. Indeed, 
a naive geometric estimate of the a 4 coefficient leads to a contribution of the order of 1000a 4 JS] , which is of the 
same order as the one of A 2 , and then justifying the arguments which motivate its introduction as a model for 
the unknown higher order terms. 
• The D = 4 contributions read [23 El El : 
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where the last term is due to the A 2 term (30], and £(3)=1.202... 
• The D = 6 contributions read UJ: 



(19) 



(20) 



where p ~ 2 — 3 parametrizes the deviation from the vacuum saturation assumption of the four-quark condensate. 
We shall use as input A 3 = (375 ± 25) MeV for three flavours and [31)1 I3T1 1371 15] : 

(m„ + m d )(uu + dd) = -2m 2 / 2 

(m s + m u ){ss + uu) ~ —2 x 0.7™^/^ 

(ss)/(uu) ~ 0.7 ±0.2 

a s A 2 ~ -(0.07 ± 0.03) GeV 2 

(a s G 2 ) ~ (0.07 ±0.01) GeV 4 



pa s \uu) 



(5.8 ± 0.9) x 10" 4 GeV 6 



(21) 



where: / ff = 93.3 MeV, fx — l-2/ ff . We have taken into account a possible violation of kaon PCAC as suggested 
by the QSSR analysis PAIS], 



4 Parametrization of the spectral function 

• For the resonances, we parametrize the spectral function within a narrow width approximation (NWA) 4 by 
using the most recent data compiled in PDG [35] for the </>(1019.7) and ^'(1680) with: 

r0(ioi9.7)^ e+e - =i (1-27 ± 0.02) keV , r 0(168O) ^ e+e - ~ (0.43 ± 0.15) keV . (22) 

• For the non-resonant contributions in the region below y/t < 1.3 GeV, we use the sum of the exclusive rates 
of the I — channel compiled in and use a SU(3) symmetry for keeping the ^-component. An analogous 
parametrization has been used successfully for the accurate estimate of the hadronic contribution to the muon 

4 A parametrization using a Breit-Wigner form leads within the errors to the same result. 



anomalous magnetic moment |41| . 

• Above the (/>(1680), we use a QCD parametrization of the spectral function as: 

ilmn ss (£ > t>) ~ 6{t - t>) -U 1 + a s (i)(l + 2m 2 Jt) + 1.6398a 2 (t) - 10.284a^(t) I , (23) 
where t> is the QCD continuum threshold. 



5 FESR 

We shall use the combination of FESR introduced recently in j20] for extracting m s from the V + A component 
of the t decay data 5 : 

Sxo = Mo- = / ° dt ( 1 - 2p\ -Imn ss (t) , (24) 

tc Jo \ t c J IT 

which is sensitive, to leading order, to m 2 and A 2 . Unlike the individual sum rules, these combinations of sum 
rules are less sensitive to the high-energy tail of the spectral functions (effect of the i c -cut), as it is chosen such 
that at t = t c the integral vanishes. 



Test of duality 

In principle, the value of the t c -cut of the FESR integrals is a free parameter. We fix its optimal value by looking 
for the region where the phenomenological and QCD sides of the ratio of moments: 

ftio = 2^1 , (25) 
Mo 

are equal. We present this analysis in Fig. ^ by showing the value of t c predicted by the sum rule versus t c 
and by comparing the result with the exact solution t c = t c expected to hold for all values of t c because the 
continuum is paramterized by the QCD expression above t c . From Fig. 1, one can deduce, that QCD duality is 
best obtained at: 

t c ~ (6.0 ± 0.5) GeV 2 . (26) 




Fi gure li FESR prediction of t c versus t c in GeV 2 . The green curve corresponds the central value of the data; the blue ones 
to the larger and smaller values of the data; the red line is the solution t c = t c , where both beyond this value, the predictions 
from the LHS and RHS should exactly coincide. The curves correspond to the value t> = 1.75 GeV 2 ; rh s = (56 — 145) MeV and 
a a \ 2 = -(0.04 - 0.1) GeV 2 . 



where one expects to get the optimal value of m s from FESR. In order to get this number, we have used the 
value of the invariant mass rh s = (56 — 145) MeV, which is like A 2 a tiny correction in this duality test analysis. 
Once we have fixed the value of t c where the best duality from the two sides of FESR has been obtained, we can 
now estimate m s . 

5 More technical details of this sum rule can be found in 1201 . 



Estimate of m s versus A 2 



In principle, t >1 the beginning of the QCD continuum parametrization is a free parameter. We study the stability 
of the m s output against the i>-variation and found that in the range \/t> = (1-75 ± 0.05) GeV, m s is quite 
stable. We present the result of the invariant mass rh s for different values of A 2 in Tabled 



Table 1: m s versus A 2 for t c = 6 GeV 2 and £> = 1.75 GeV 2 . 



-a s A 2 in GeV 2 m s in MeV 



0.02 77 ±55 ±15 ±8 

0.04 94 ±36 ±10 ±7 

0.06 108±30±9±6 

0.07 114±27±9±6 

0.08 120±25±6±5 

0.10 131±22±5±4 

0.12 142±21±5±4 



The 1st error is due to the data, the 2nd one to t c and t> , and the 3rd one to the QCD parameters. 

Using the value of a s A 2 given in Eq. 1)2 H|. we deduce the predictions: 

m s = (114±27±9±6± 19) MeV => m s (2 GeV) = (100 ± 28 cxp ± 20 t h) MeV , (27) 

where the last error in rh s is due to A 2 6 . We have used the conversion scale: 

m s (2 GeV) ~ 0.876 rh s . (28) 

One can notice that the result is perfectly consistent with the one from r-decay (201 given in Table |31 using the 
same combination of FESR. Like in the r-decay data, the theory with A 2 = tends to give too small value of m s 
though the error is quite large. 



6 r-like sum rules revisited 

r-like sum rules have been proposed in SN |15| for extracting m s from e + e~ data, which have been exploited 
later on in |16U17| using the same or a slight variant of the SN-sum rule. SN JS] results have been confirmed by 
[T7] using the expected small effects of the SU(2) breaking due to u> — p mixing. However, the central value of 
the results obtained in ^1 El are slightly higher than recent estimates, though consistent within the errors. In 
the following, we reconsider the original sum rules: 

^'* = ^r Sjsw l ds ( 1- ^) + w CTe+e -^' -' (29) 

and the S'C/(3)-breaking combinations: 

Ai^ = R T .\ — Rr,4>- (30) 

Here, Sew = 1.0194 is the electroweak correction g^l and |Kd| 2 = 0.975 is the CKM mixing angle. The QCD 
expressions of these sum rules have been given in |15j . except that we shall replace the contributions of the 
uncalculated of and HO terms of the PT series by the contribution of A 2 . We shall also use the computed 
coefficient of the aj?m 2 term — 250 .4 [2E1 instead of the estimated 218.55 used in Contrary to the 

previous sum rules, the r-like sum rule is more precise near the real axis due to the presence of the threshold 
factor (1 - s/A/ 2 ) 2 . 

6 We have taken the central value of the asymmetrical errors due to A 2 . 



Upper bound on m s from R T ^ 



We shall use the positivity of R T ,<f> and saturate the spectral function by the 0(1019.7) contribution. In this way, 
we derive a lower bound on R Ty ^, which we show in Fig. Inversus M T . 




\ 3 



l .2 
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Figure 2: Upperbound of Rr <f> versus M T in GeV for a s X 2 = —0.07 GeV 2 , using the central value of the data. 

From this figure, one can derive: 

m s < (172 ± 16) MeV => m s {2 GeV) < (151 ± 14) MeV . (31) 

This upper bound is comparable with the previous value m s (2 GeV) < (147 ± 21) MeV obtained in 15 using 
the same method, where the improved accuracy comes mainly from a more precise value of the r^ e + e - width. 
This bound is also comparable with the one 148 MeV from a direct estimate of the quark chiral condensate |26| . 

Extraction of m s from R T ^ 

Parametrizing the spectral function by the 0(1019.7), 0(1680) and the non-resonant contributions below 1.39 
GeV, we deduce, for the QCD parameters in Eq. J2U, the results in Table [3 Comparing these results with 
the previous ones in |15|. we notice that the inclusion of the non-resonant states contributions has increased 
slightly the value of i? T) 0. However, this small change affects the value of m s which is also a correction in 
the QCD expression of R r .$. Taking as an optimal estimate of m s the value which is stable in the change of 

Table 2: Phenomenological estimates of R T j and central values of rh s to order 



1.4 2.82 ±0.47 

1.6 1.73 ±0.05 111 

1.7 1.66 ±0.05 128 

1.8 1.65 ±0.06 129 

1.9 1.65 ±0.06 127 
2.0 1.80 ± 0.11 



M T ~ (1.8 ± 0.1) GeV 7 , we obtain: 

rh s = (129 ± 15 oxp ± 25th) MeV =^ m s (2 GeV) = (113 ± 13 oxp ± 22 th ) MeV . (32) 

7 Notice that the inclusion of higher states has slightly shifted the position of the optimum from 1.6 GeV (Fig. |5J to 1.8 GeV 
(Table |5J . 



This value is consistent with the result in Eq. (|27J) . Like the previous FESR, this prediction is also affected to 
leading order by A 2 , where the value of m s increases with — a s X 2 . For A 2 = 0, the corresponding value of rh s is 
104 MeV, which is still consistent with the other determinations, though in the lower side. Therefore, one can 
notice that, contrary to the r-sum rule, this sum rule cannot differentiate between the two cases A 2 = and A 2 7^ 0. 



Extraction of m s from Ai^ 



Here, we analyze the Ai^ sum rule. Unlike R T ,<p, ^1,4, is not sensitive to leading order to A 2 . However, like the 
S'[/(3)-breaking sum rule, it has the inconvenience to involve the difference of two large independent channels 
(isovector-isoscalar here). Using the optimal value of R T .4> at M T = 1.8 GeV from Tabled and using the value 
of R T .\ — 1.78 ± 0.025 at this scale from r-decay data ^SIi one can deduce: 

A10 (1.8 GeV) =0.13 ±0.06 , (33) 

leading to: 

m s = (113 ± 26 oxp ± 4th) MeV m s {2 GeV) = (99 ± 23 oxp ± 4 th ) MeV , (34) 

which is in good agreement with the former determinations. 



Final value of m s from e + e 



One can see in Table that there are good agreement between results from e + e~ data from alternative forms 
of the sum rules, indicating the reliability of the results. Taking the average of the results in Eqs. (|27|l . I|32(l and 
()34p. we deduce the final value from e + e~ data: 

rh s = (119 ± 17) MeV =^ m s (2 GeV) = (104 ± 15) MeV . (35) 

This value is consistent within the errors (though in its higher side) with the old value: m s (2 GeV) =(125 ± 
14 oxp ± 20th) MeV in [TJ], and the one (139 ± 31) MeV in [T7]. 



7 Present status of light quark masses 



In this part, we present the status of the recent determinations of the light quark masses from different sum 
rule channels, to the same order and including the A 2 term, which we also compare with the lattice results 
including dynamical fermions. 



The (pseudo) scalar channels 

These channels are, in principle, the best place for extracting the value of the light quark masses because these 
masses enter as the leading overall coefficients in the corresponding QCD correlators, precisely known up to order 
ol\ or alternatively to order o? s plus the A 2 -term which mimics the unknown higher order terms. 
• Estimates of the sum of the light quark masses (m u + wid) [8] have been updated in |3U1 |S] by including 
the A 2 and ai? corrections, and by using the parametrization of the 37r spectral function |S] which satisfies the 
ChPT constraints. The result is: 



(m u + m d )(2 GeV) ~ (8.6 ± 2.1) MeV , (36) 
where A 2 decreases the value of the sum by 5%. Combined with the ChPT ratio 0: 

r 3 = 2 ™ s =24.4 ±1.5, (37) 
(m u + m d ) 

one can deduce the value of m s given in Tabled 8 - 

• Direct extractions of m s from kaon sum rules also exist in the literature Here, the analysis suffers 

from the unmeasured value of the kaon radial excitations -ftT(1460) and -ftT(1830) decay constants which play an 

8 Kaon sum rule also gives an analogous value UHHE1 but the corresponding spectral function is less controlled by ChPT than 
the one of the pion. 



important role at the scale where the sum rules are optimized. We expect that the errors induced by this model 
dependence have not yet been properly included in the quoted small errors of the estimated decay constants. 

• The scalar channel has been revisited in ^2] using Ktt phase shift data 9 , with the resulting value of the quark 
mass given in Table [3] Here, the phenomenological side of scalar sum rule is better known due to the availability 
of the i^7r-phase shift data combined with the constraints from ChPT. 

• Lower bounds on the light quark masses have been also derived from the (pseudo)scalar channels using the 
positivity of the spectral function [25] . These bounds have been updated in [5] by including the effect of A 2 
and the order PT contributions. The best updated bounds from pseudoscalar kaon and pion sum rules given 
in including the and l/g 2 -term are: 

m s (2 GeV) > (71 ±4) MeV , {m u + m d ){2 GeV) > (5.9 ± 0.3) MeV . (38) 

The inclusion of the ol\ term obtained in 03] decreases this value by about 2 MeV. However, slightly different 
central values without error bars have been given in |44j but agree with ours within our quoted errors. 

r-decay data 

r-decay data have been used using different methods for extracting m s : 

• In |181 119j . S'C/(3)-breaking moment sum rules involving the difference of the non-strange and strange V + A 
components of T-decay data have been used. These sum rules have the advantage to be non-affected by A 2 to 
leading order, but have the inconvenience to involve a strong cancellation of the two independent channels ud 
and us. The value of m s is given in Table |2| 

• Alternatively, a combination of FESR involving only the AS = — 1, but sensitive to A 2 to leading order, has 
been proposed in |20|. This sum rule has been used for studying the effect of A 2 on the value of m s . Giving the 
value of A 2 in Eq. I|21|l , one obtains the value in Table 01 

9 For some recent discussions on the scalar channels from the sum rules, see e.g. 1451 . 



Table 3: Recent phenomenological determinations of m s (2 GeV) to order a 3 including the tachyonic gluon mass 
A 2 which parametrizes the UV renormalon contributions into the PT series. 



Channels 


Refs. 


m s (2 GeV) in MeV 


e + e data 

FESR 

Rt,4> 

Average 


This work 
This work 
This work 
This work 


100 ± 28 cxp ± 20 th 
113 ±13 cxp ± 22th 
99 ± 23 oxp ± 4th 
104.3 ± 15.4 


r-decay data 

FESR 

SU(3) breaking SR 


m 

[HO EH 


93 ±30 
81 ±22 


Pseudoscalar 

l lull 0-t\ v_/iiir .1 

Kaon FESR 

Kaon Exponential SR 


jom in |oj 

D 

CDI+) 


100 ±12 
103 ±9 


Scalar 

K* SR 


[12 * 5 


88 ±8 


Chiral condensate 

N , B* -B,D^ Klv 




131 ±18 


Final Average <> 


Weighted 
Arithmetic 


96.10 ±4.80 
96.30 ± 17.5 



Result at order of quoted in 1101 : of corrections increase the value by 2 MeV. 
*> Result to order of quoted in 1121 . of correction is expected like in the pseudoscalar channel to be small. 
**) Not included in the average as known to order a s . 

t) We have assumed that the different determinations are independent from each others. 



Direct extraction of the chiral condensate 

Extraction of the chiral condensate (ipip) nas been used in |26| for estimating and bounding m s . The lower and 
upper bounds The nucleon and B* — B sum rules give |26| : 

(ipip)(M N ) ~ [-(225 ± 9) MeV] 3 , (39) 

which combined with the GMOR relation and the ChPT mass ratio in Eq. (|37|l leads to: 

{m u + m d )(2 GeV) ~ (10.8 ± 1.3) MeV, and m s {2 GeV) ~ (131 ± 18) MeV . (40) 

Bounds have been also derived from the D — > K*lv decays |26| : 

0.6 < (^V)(l GeV)/[-229 MeV] 3 < 1.5 , (41) 

giving: 

6.8 MeV < (m u +m d )(2 GeV) < 11.4 MeV . (42) 
Combined with the ChPT mass ratio r^ in Eq. (|37|) . it gives: 

82 MeV < m s {2 GeV) < 138 MeV , (43) 



which are comparable with the lower bound from the 
e + e~ data given in Eq. I|31fl . However, as the result is 
our final estimate. Instead, we use the allowed region 



pseudoscalar SR given in Eq. I|38|l and upper bound from 
obtained to order a s , we will not consider these bounds in 
in order to deduce the inaccurate estimate given in Table 



Final value of the strange quark mass and m b /m s from QSSR 

• As a final result, we consider the weighted average given in Table which emphasizes the contributions of the 
most accurate results from (pseudo)scalar channels, which give more weight in the averaging procedure: 

m s {2 GeV) = (96.10 ± 4.80) MeV . (44) 

• As discussed in previous subsection, the precisions from these two channels can be qualitatively understood 
because the square of the strange quark mass enters as the leading overall coefficient in the analysis of the 
correlator associated to the divergence of the axial-vector [resp. of the vector] currents, while in the vector and 
tau-decay channels m s enters as ml/q 2 corrections in the corresponding two-point correlator. The accuracy of 
the phenomenological side of the pseudoscalar sum rule is more questionable due to the lack of data and to the 
accuracy of the radial excitation decay constants which play a crucial role in the analysis. The phenomenological 
side of scalar sum rule is in a better shape due to the availability of the K-K--pha.se shift data and to the constraints 
from ChPT. A confirmation of the accuracy obtained from the (pseudo) scalar channels requires an independent 
analysis of these channels. 

• However, it is difficult to quantify with a good precision the systematic errors of the different sum rules 
approach, though, in each analysis, the different authors have used their own estimate of such errors by studying 
the effects of external parameters (sum rule scale, continuum threshold....) based on optimization or/and stability 
procedures or duality tests. Due to the remarkable good agreement of the different results given in Table|3within 
about 1 cr, we might expect that the quoted error in Eq. Q44JI is quite realistic though relatively small. This value 
of m s , is consistent with the older sum rule (arithmetic) average (117.4 ± 23.4) MeV quoted in though in 
the lower side, indicating the stability of the sum rule results with time and then their reliability. 

The more conservative value from the arithmetic average obtained in Table |3| 

m s (2 GeV) = (96.3 ± 17.5) MeV, (45) 

can be translated into the range of m s values allowed by the sum rules analysis : 

79 MeV < m s (2 GeV) < 114 MeV, (46) 

which can be compared with the rigorous lower [resp. upper] bound coming from the positivity of the spectral 
functions in the pseudoscalar [resp. <fi] sum rules updated to order 10 : 

(71 ± 4) MeV < m s {2 GeV) < (151 ± 14) MeV. (47) 

• The final average value of m s obtained from phenomenological methods given in Eq. (|44|l are inside the range 
of values quoted by PDG |39| : 

m s (2 GeV) = (80 - 130) MeV . (48) 

• One can also compare the value in Eq. (|44|1 with the different lattice results [^I^UEIl- The recent lattice 
results including dynamical fermions are: 

m s (2 GeV)| n/=2 = (100 ~ 130) MeV , m s {2 GeV)| n/=2 +i = (70 ~ 90) MeV , (49) 

which appear to depend on the number of flavours. The difference of the results for rif = 2 |22| and nt = 2 + 1 
|23"| (see however j^J) and the slightly higher prediction of the ChPT mass ratio r 3 = 27.4 ± 4.2 defined in Eq. 
(|37|l may indicate that, it is premature, at present, to extract a precise value of m s from the lattice calculations, 
before a reliable control of the systematic errors, higher order terms and some other effects. 

• Running m s until m^mi,) = (4.23 ± 0.06) GeV [351|5J^211i by taking care on the threshold effects, one can 
deduce the useful scale-independent quantity for model-buildings: 

r 5 = — = 50 ± 3 . (50) 



10 Stronger lower and upper bounds from a direct extraction of the chiral condensate has been obtained in Eq. 1431 . but they are 
only known to order a s . 



Implied values of the up and down quark masses from QSSR +ChPT 



• Using the previous value of m s (2 GeV) in Eq. fZty together with the ChPT mass ratio in Eq. 1(37(1 . we can 
deduce: 

(m u + M d ){2 GeV) = (7.9 ± 0.6) MeV , (51) 

in nice agreement with the result from the pion sum rule in Eq. 1(36(1 . Using again the ChPT mass ratio 11 : 

ill 

— = 0.553 ±0.043 , (52) 
m d 

we obtain: 

m d {2 GeV) = (5.1 ± 0.4) MeV , m„(2 GeV) = (2.8 ± 0.2) MeV . (53) 

• Taking the average value of m s (2 GeV)= (95.8 ± 4.9) MeV, by excluding the pion sum rule result in Table 
and using the prediction of (m u + md) from the pion sum rule in Eq. 1(36(1 ■ one can deduce the ratio: 

2t7£ 

r 3 =- -2-^ = 23.5 ±5.8, (54) 

(m u + m d ) 

in perfect agreement with the ChPT mass ratio in Eq. 13 7|) . 



8 Conclusions 



We have revisited the estimate of the strange quark mass from e + e~ data. Including the plus a phenomeno- 
logical estimate of the UV renormalon contributions parametrized by the tachyonic gluon mass A 2 , we deduce 
the final value from e + e~ data in Eq. 1(35(1 and reported in Tablc|3J We compare this value with recent determi- 
nations from different channels in Table known to the same level of approximations. Our final result coming 
from a weighted average of different determinations from Table |3] is given in Eq. (|44l) . The updated lower and 
upper bounds for the strange quark mass to order are summarized in Eq. (|47l) . The value and range of m s 
given in Eq. H44|) to Eq. (|47|) are inside the PDG values quoted in Eq. I|48l) and agree within the errors with the 
recent lattice calculations in Eq. (149(1 including dynamical fermions. Combining the final average result of m s in 
Eq. (|44|l with the ChPT mass ratios, we deduce the value of the running u and d quark masses in Eqs. 1(51(1 and 
lj53"|l . while we also predict in Eqs. (|54*jl and l(5T)| . the useful scale-independent mass ratios m s /(m u + md) and 
m,b/m s . 
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